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Abstract The phase-space structure of our Galaxy holds the key to understand and reconstruct its formation. 
The ylCDM model predicts a richly structured phase-space distribution of dark matter and (halo) stars, 
consisting of streams of particles torn from their progenitors during the process of hierarchical merging. 
While such streams quickly loose their spatial coherence in the process of phase mixing, the individual stars 
keep their common origin imprinted into their kinematic and chemical properties, allowing the recovery of 
the Galaxy's individual "building blocks". The field of Galactic Archeology has witnessed a dramatic boost 
over the last decade, thanks to the increasing quality and size of available data sets. This is especially true for 
the solar neighborhood, a volume of 1-2 kpc around the sun, where large scale surveys like SDSS/SEGUE 
continue to reveal the full 6D phase-space information of thousands of halo stars. In this review, I summarize 
the discoveries of stellar halo streams made so far and give a theoretical overview over the search strategies 
imployed. This paper is intended as an introduction to researchers new to field, but also as a reference 
illustrating the achievements made so far. I conclude that disentangling the individual fragments from which 
the Milky Way was built requires more precise data that will ultimately be dehvered by the Gaia mission. 

Keywords Milky Way: dynamics — Milky Way: kinematics — Milky Way: solar neighborhood 



1 Introduction 

Not until the past two decades have we begun to understand the composition of our Galaxy in a universal 
context of galaxy formation and evolution, thanks to great progress in both observations and theoretical work. 
Historically, there have been two differe nt models for the early formation of the Milky Way's stellar halo. 
Eggen. Lvnden-Bell and Sandage I (119621) found strong correlations between the metallicity and eccentricity 
of stars in the Solar neighborhood. To explain the highly eccentric orbits of the oldest stars, they argued that 
at the time of formation of these stars, the Galaxy could not have been in dynamical equilibrium, but in a 
state of rapid gravitational collapse (~ 10^ yr) from a larger homogeneous spheroid, the "protogalaxy". After 
the Galaxy reached dynamical equilibrium and became rotationally supported, further star formation would 
have taken place in a metal-enriched disk, thereby explaining the disk-like orbits of the metal-rich stars. This 
sc enario is called a monolithic collapse. Alth ough the correlation between eccentricity and metallicity found 
bv lEggen. Lvnden-Bell and Sandage I d 1962b is empirically correct for samples including both disk and halo 
stars, many authors later questioned its validity for a pure h alo sample and it was de finitively discarded from 
a kinematica lly unbiased samp le of 1203 metal-poor stars ( IChiba and Beers Il2000l and references therein). 
In 1978, ISearle and Zinn I published their seminal study on the distribution of globular clusters in the 
outer Milky Way halo; they found no abundance gradient, but a substantial spread in the ages of the clusters. 
They concluded that, although the lack of an abundance gradient would still be consistent with cluster forma- 
tion during the rapid free-falling phase of a galactic collapse (because the freefall timescale is much shorter 
than the enrichment timescale), globular clusters would not have su ch a high internal chemi cal homogeneity, 
if they would have formed during a monolithic collapse. Therefore. ISearle and Zinn I (Il978h proposed a slow 
(<: 10^ yr) and more chaotic process in which the Milky Way was gradually built through the merging of 
several small protogalactic "fragments". This is called a bottom-up scenario, in which galaxi es are generally 
formed through the amalgamation of smaller fragments. However. IChiba and Beers I (120001) found that the 
density distribution of the halo changes from highly flattened in the inner parts to nearly spherical in the 
outer parts with no discrete boundary. Also they found a continuous negative rotational velocity gradient 



with height above the Galactic plane, both results that are hard to arrange with the original [SearleandZinn 



(11978.) scenario, because chaotic merging would not produce such kinematic and spatial structures 



Therefore both simplistic models cannot fully explain the distribution, kinematics and abundances of 
halo stars in the Milky Way. However, the stellar halo only accounts for ~ 1 % of all stars in the Milky Way, 
and a more general approach towards understanding the formation of the Milky Way, and disk galaxies in 
general, is needed. Such an approach has to take into account cosmological principles and the nature and 
properties of dark matter. Dark matter is the dominant form o f matter in the Universe and it should play 
the ke y role in galaxy formation through gravitational collapse ( IPress and Schechterlll974tlWhite and Rees 
19781) . Numerical and semi-analytic models show that, if the random velocities of the dark matter particles 
are small compared to the speed of light (cold dark matter or CDM), present day galaxies can be formed 
through gravitational collapse of the small density fluctuations which where present in the early Universe 
after the big bang. While the average density of the cosmic matter declined as the Universe expanded, 
some density enhancements of sufficient size became more pronounced, because of their own gravitational 
attraction. These protosystems drew in more matter from surrounding regions and subsequently merged to 
larger and larger structures, increasing further the lumpiness of the once highly (but not perfectly) uniform 
distribution of matter. The "trunks" of such "merger trees" in the hierarchical scenario are the present day 



galax ies (e.g.'Katz lll992l: ICole et al I1994|:IMo et al. ll998HSteinmetz and Navarroll999l:lBekki and Chiba 
200 U Steinmetz and Nav arro Il2002h . Although these simulations lack high enough resolution to study the 
detailed spatial and kinematic distribution of stars, they qualitatively can explain the formation and properties 
of different galactic components. Recent simulations of galaxy formation in the ACD M cosmology were ab le 
to build Milky Way-like spiral galaxies including thin disk , thick disk, bulge and halo JAbadi et al. Il2003allbt 
Brook et al. ll20o4lSamland1l2004l;lGovernato et al.1l2007h . 

The hierarchic al merging scenario pr edicts a richly structured phase-space distribution of dark mat- 
ter and stars (e.g. iBullock and Johnston 2005b. Early e vidence for remnants of accretion events in the 
outer halo was found by, e. g., iDoinidis and Beers I (Il989[ five sta r s at (d ) ~ 5 kpc distance from the sun), 
Arnold and Gilmore I (Il992[ four stars at t/ ~ 30 kpc), llbata et al.l ( 11994 discovery of the disrupting Sagit 



tarius dwarf galaxy as a clump in radial velocity space of several 100 K and M giants) or lMaiewskiet al. 



(Il99d several metallicity/ phase-space clumps in a sample of 154 stars at li < 8 kpc). The last years have 
seen a dramatic increase in the detection of halo substructur e thanks to automatic large-scale photomet- 
ric surveys such as 2MASS JSkrutskie et al. 1120061) and SDSS JStoughton et al.ll2002l) . Prominen t examples 



are the currently accreted and disrupted Sagittarius dwarf galaxy with it s associated tidal tails (llbata et al. 



200ll:lMaiewski et albOOllBelokurov et al. l2006al). t he Orphan stream (jBelokurov et al. Il2006al) . die Virgo 

overdensitv dVivas et al. I 200ll: iNewberg et al. 112002c Ijuric et alj|2008h or the globular clusters Palomar 5 

dOdenkirchen et al. l200ll:lGrillmair and Dionatos l2006h and NGC 5466 ( JBelokurov et al. l2006bl:lGriUmair and Johnson 
2006h . 

Numerical simulations predict that mo st accreted satelhtes wou ld also spread their tidal deb ris on eccen- 
tric orbits into the Milky W ay's inner halo (IHelmi and White Il999h or even the disk component dAbadi et al. 



2003bclMeza et al. l2005b : thus we can expect to find such "fossil remnants" right in our immediate surround- 
ings. Although no longer spatially coherent, such stellar streams keep their common origin imprinted into 
their chemical and dynamical properties. This allows us to study the formation history of our Galaxy by 
disentangling its individual "building blocks" from the phase-space distribution of solar neighborhood halo 
stars, for which we have much better kinematic and abundance informations than for their distant counter- 
parts. The difficulty, however, lies in the very low density contrast that is expected for individual streams. 
Due to the short dynamical time scales near the disk, a single progenitor can dispose multiple streams in 
the solar neighborhood in a sho rt period of time, leading to an expected total of several h undreds of cold 



tne solar neignbornooa m a sno rt penoa oi tune, leaumg to an expecteu total oi several n unareas or com 
(o- < 5 km s ') steflar streams (iHelmi and White lll999riHelmi etaP 120031: lGouldll2003lPl . For obtaining 



' C ompare this number to an expected number of tens of thousands of CDM particle streams JHelmi et al. l2003l ; IVogelsberger et aJTI 

Eooi). 



most of the informations hidden in the phase-space structure of local stars it follows that (i) a large sample 
of thousands of halo stars with very accurate velocities would be needed and (ii) the strategies for finding the 
individual streams and deriving informations about the past should be optimized. Th e first point is currently 



addressed by large-scale surveys such as RAVE (Steinmetz et al. 2006) and SEGUE dYannv et al. P 2009) and 
will be completely satisfied when the Gaia mission jPerrvman 11200 lb is completed. The second point is part 
of this review, which should be seen both as an introduction for people new to this field as well as a reference 
summarizing the current state of research. 



2 The formation of stellar halo streams 

In this section I will summarize the main processes that lead to the formation of stellar halo streams that 
might show up in samples of nearby halo stars. 



2.1 Tides and dynamical friction 

Two effects are important for accreted satellites to reach the inner halo and the disk. The first is dynamical 
friction. As the satellite moves through the fluid of lighter particles of the host halo, they get pulled towards 
the satellite, which in the meantime moves along. The resulting overdensity or "gravitational wake" behind 
the satellite leads to a steady deceleration of the latter. The deceleration of a satellite of (point) mass Ms 
moving at velocity v, on a straight orbit through an infinite, homogeneous and iso tropic field of particle s 
with mass density phost can be expressed through Chandrasekhar's analytic formula (IChandrasekhar II 1 943h : 



^ = -4nG^M,lnApi,U< Vs)-^. (1) 

at I'UsI'' 

Here, G is the gravitational constant, pho,st(< Vs) the density of field particles at the position of the satellite 
with velocities less than \Vs\ and In A is the Coulomb logarithm, named so in analogy with the equivalent 
logarithm in the field of plasma physics. Note that the drag force, Ms-jf-, is proportional to M^, so that 
massive satellites loose their energy and angular momentum much faster than low massive ones. Equa- 
tion [H further predicts the m erging timescale to depend on the masses of the satellite and host halo as (e.g. 



Binney and Tremaine OOOSl eq. 8.13) 



Tmerge oc (inA— -^) . (2) 

Besides assuming linear motion of a point mass through an uniform and isotropic background of particles, 
Chandrasekhar's formula also neglects any gravitational attraction of these particles for one another There- 
fore, in general, equations ([Tji and (|2|i can not describe accurately the orbital decay of extended, mass-loosing 
satellites moving on eccentric orbits within a live host halo in A^-body simulations. For example, while 



Sales et al. I (l2007h found that the infall timescales of satellites onto host halos in their A^-body/gasdynamical 
simulations follow the trend with mass ratio given by equation (|2]i, Boylan-Kolchin et al. (2008) found a 
non-linear dependence in their A^-body simulations scaling as (Mj/Mhost) ''^ (i-C- the satellite sinking more 
slowly than predicted by the standard formula). Another uncertainty is the choice of the Coulomb logarithm. 
The term In A is given by 

InA-ln(^), (3) 



where /7min and /^max denote the maximum and minimum impact parameters, i.e. the maximum and minimum 
distances of closest approach between the satehite galaxy an d the field star s of the host galaxy. Usually, bmin 
is set to the size of the satellite ga laxy (i.e. its virial radius. White ^4976 j), while the choi ce of /j^ax is less 
clear and somewhat arbitrary (Binn ev and Tremaine 1 12008, § 8.1). Hashi moto et al. I (l2003b argued that /j^ax 
should be varied according to the actual separation of the satellite from the host g alaxy's center in order to 
account for the density gradient of the host halo in A^-body simulations. However, Fuiii et al. I (l2006h showed 
that this approach predicts too fast sinking times of satellites which suffer significant mass loss, because 
recently stripped particles would contribute to the drag on the satellite. 

Mass loss of the satellite is the second important effect and occurs primarily at each pericentric passage 
through tidal shocks, but also at other parts of the orbit due to the limitation of bound satellite particles to stay 
within a tidal surface. If the satellite is on a circular orbit, Jacobi-Hill theory predicts that stars with Jacobi 
energy Ej - jv^ + ^eff stay bound to the satellite, if they lie inside the last closed zero velocity surface. 



the so-called Roche surface, at which Ej equals the effective potential (Figure [Tl and iBinnev and Tremaine 



l2008i § 8.3). The latter is given through 0eff(r) - <l>(r) - ^Iw x rp, where 0(r) is the gravitational potential 
and the second term describes the centrifugal force with the angular speed oj of each mass around the 
common center of mass. The Roche surface includes the two Lagrange points Li and L2, which are saddle 
points of the effective potent ial and are invers ion symmetric with respect to the satellite's center only for 
very small satellite masses tChoiet al.ll2007h . Stars bound to the satellite and having Jacobi energy Ej 



can escape the satellite, if weak perturbations elevate Ej above 0eff(Li) or 0eff(L2). So, mass loss occurs 
near the two Lagrange points, leading to the formation of a leading and trailing tidal t ail. Such tails can 
be observed as a spectacular example stretchi ng from the Sagittarius dwarf spheroidal flb ata et al 



Majewski et al. 2003; Beloku rov et al. Il2006ah or the globular cluster Palomar 5 (Odenkirch en et al. 



2001 



2001 



Grillmair and Dionatos II2006I) . The leading tail contains stars at lower energies than the trailing tail, so each 



time stars get stripped from th e satelhte, e.g. through tidal shocks at pericenter, roughly two gro upings of 
stai-s at distinct energies resuh (|johnstonlll998HMeza et al. II2005I; IChoi et al.1l2007b . l¥>^rnick et a l. (200^ 
proposed that this energy difference could be used for the distinction b etween lead i ng an d trailing stars 
even after significant dispersion throughout the Galaxy. As pointed out bv lChoi et al. I (120071) . (dark matter) 
particles in the leading tail continue to be decelerated by their satellite's gravity and fall towards the central 
regions of the host. It has not been studied yet, however, if this could provide another mechanism for stars 
to populate the inner halo, even if their satellite stays hung up in the outer halo. 



2.2 The contribution to the inner halo 



Both the adiabatic growth of the host halo and the energy loss of the satellite through dynamical friction have 
the effect of its orbit spiraling in towards the Galactic center. There is still debate whether th e orbits also 



circularize during infa ll, and if so, whether dynamical friction causes this circularization (e.g. IColpi et al. 



1999t IgIU et al. |2004|) . However, the decay of the satelUte's orbit leads to a stronger tidal field and in c reases 
the mass loss, which in turn decreases dynamical friction. Using a fully analytical model, IZhao I ( 120041) 



estimated this interaction between dynamical friction and mass loss to delay the merging of typical present 
day satelhte s (10^ - lO^M©) with initial distances d > 20 kpc to more than a Hubble time. In support of this. 



Colpi et al. I (119991) already argued that the effect of mass loss outperforms dynamical friction and makes it 
unlikely that remnants of satellites with Mj/Mhost > 0.1 have reached the inner halo yet. 

Numerical simulations of the hierarchical buildup of Milky Way-like galaxies in the ACDM context 
seem to support these dynamical friction arguments. A consistent finding is that the inner halo consists of 
particles stripped from a few significant contributors. These are more massive and have been accreted much 




Fig. 1 Contours of constant effective potential in the restricted three-body problem (Jacobi-Hill theory) 
for satellite-to-host halo mass ratios m/M = 1/5. Note the non-circular form of the Roche surface, which 
contains Li, and the asymmetry of Li and L2 with respect to m. Stars with Ej < 0eff(Li) are bound to the 
satellite. 



earlier (>8 Gyr ago) than the satellites that remain today ( Bullock and Johnston I l2005t Font et al. I l2006t 
Sales et al. II2007I : iKazantzidis et al. Il2008t iRead et al. II2OO 8), although the possibility exists that disruptio n 
of the massive satellites is not yet completed and remnants still exist in the inner halo (.Cooper et al. Il2010r) . 
Stated otherwise, "the 'building blocks ' of the stellar halo were on average more massive and w ere accreted 



and d isrupted much earlier than the population of satellites that survive until the present." (ISales et al. 



20071) . So, in a cosmological context, the presently observed satellites are not assumed to be representative 
for the ones already accreted into the inner halo. Therefore it is not surprising that present-day satellites 
exhibit different che mical compositions than halo stars. The latter tend to have higher [a/Fe] ratios at a given 
[Fe/H] (iVenn et al. 12004 ). where a stands for the average abundance of the elements Mg, Si, Ca and Ti. This 
points towards rapid star formation in their progenitor systems, because the a elements are mainly produced 
in type II supernova (SN) explosions on a short timescale (~ 10^ yr), while only after a longer timescale of 
~ 10^ yr SNela start to lowe r [a/Fe] again through expelling additional Fe into the ISM. Recently, however, 
Nissen and Schuster I ( 120101) reported a dichotomy of nearby halo stars {{d) =115 pc) with respect to [c/Fe], 
with the high-Q- stars being more gravitationally bound than the low-ff ones. This could indicate that low 
mass dwarf galaxies have contributed to the inner halo, although to what extend is still to be determined. A 



possible mech anism for dwarf galax ies to reach th e inner halo is infa ll in massive, loosely bound groups, as 
pointed out bv lLi and Helmi I (l2008b . According to lRead et al. (120081) . this allows even low mass satellites to 
merge into the inner halo in less than a Hubble timeJ Read et al. I also have shown that - assuming isotropic 
infall - a typical MiUcy Way halo will have one out of three mergers at a given mass within 20° inclination to 
the disk plane. While the high -inclination satellites would primarily contribute to the Milky Way stellar halo 
dBuUock and Johnston 112005b . the low-inclination ones would be dragged completely into the disk plane 
by dyna mical friction (Fig. 6 in Read et al. 2008), which results in a richly structured accreted thick disk 
(see also lAbadi et al. Il2003bt iMeza et al. II2005I) . In all cases the encounters would heat the pre-existing thin 
disloand produce several morphol ogical features like bars, flares and warps (lAbadi et al. l2003btlRead et al. 



12008; Vill alobos and Helmi Il2009|). or even overdensit ies of stars at distinct rotational velocities that travel 



as waves through velocity space (IMinchev et al. 1120091) . 

It is interesting to ask how one can distinguish between accreted and in situ thick disk stars in a metal- 
poor local sample. This would give us clues about the fraction of accreted stars versus those that were 
born within the Milky Way's potential well and contributed to the thick disk through other processes. 
Mai nly three such processes have been discussed in the literature: (i) heating of a pre-existing thin disk 
(e.g. I Villalobos and Helmi Il2008t iKazantzidis et al. 112008.) . (ii) radial mixing, i.e. migration of stars both 



from inner and outer radii th rough the solar neighborhood JSellwood and Binnev Il2002t iHavwood 



Schonrich and Binnev 1120091) . or (iii) in situ star formation during a gas rich merger (Bro ok et al 



120081 

2004. 



By directl y comparing the o rbital properties of stars from four simulations of these thick disk formation 



Sales et al. I (120091) found clear differences in the eccentricity distribution between the generated 



scenarios, 

thick disks. While accreted stars always contribute to the high-eccentricity tail of the distribution, in situ 
stars have low eccentricities e ~ 0.2 - 0.3, regardless which of the mechanisms (i) to (iii) gave them thick 
disk properties. Given t hat a certain fraction of stars in the Solar neighborhood has highly eccentric orbits 
dChiba and Beers l2000l) and assuming that they have indeed been accreted, we can further ask if we can trace 
back any of theses stars to a common origin. In other wo rds, do any of these stars belong to the same stellar 
streams? Assuming a halo solely built through accretion. iHelmi and White I ( 119991) predicted the number of 
individual streams passing the solar neighborhood to be 



A^streams~[300-500](r/10Gyr)3 



(4) 



where t is the time since disruption. In reference to the dynamical friction arguments made above, these 
streams are expected to stem mainly from a few progenitors only. The Milky Way seems to have a quiet 
merger history since f ~ 10 - 12 Gyr (corresponding to the age of the thick disk). In such a time any 
debris from an accreted satellite is expected to be spread over several tens of kpc. Therefore one cannot 
expect to find many of the halo stream s predicted by equation (|4| as spatially coherent structures in the solar 



neighborhood ( Seabroke et al.J 2008h. although hin ts for a spatial clumping of one kinematically identified 
substructure have been found (IKlement et al. 1120091) . Instead, other search strategies for stellar streams have 
to be applied, which I describe in the following section. 



^ Thereby only very massive high-redshift satellites would be able to form a thick disk from a thin disk through heating 
<Kazantzidis et al. 12003) . 



3 Search strategies for halo streams 

3.1 Integrals of Motion and Action- Angle Variables 
The collisionless Boltzmann equation 



d df ^, df d<l>df. 

_^(,,,,0=^.g(v,.^-^^) = (5) 

is the fundamental equation of stellar dynamics. It assumes that encounters between stars can be ignored 
so that every star moves under a mean gravitational potential <P(r, t) (therefore the term "collisionless"), 
where r - (x\,X2,x-i), and v - {vi,V2,v-i) denote the position and velocity vector of a star with respect 
to the center. f{r, v, t) is the distribution function describing the density of stars in phase-space similar to 
the distribution of atoms in a gas. Equation (|5]l states that the phase-space density / around the phase-space 
point of any particular star is constant. It follows that, as a sample of stars which is initially concentrated 
into a small phase-space volume (like in a satellite galaxy) disperses in configuration space, at any particular 
point in p hase-space they w i ll grow together kinematically in an infinitesimal volume around that point. For 



example. iHelmi and White I (119991) have shown that the velocity dispersions at each point along the track of 
a single stream of stars lost at the same passage decrease overall with f"' with some periodic oscillations 
due to the spatial density enhancements at the turning points of the orbit. An observer, however, will always 
measure the coarse-grained phase-space density (/), which is the average phase-space density in his finite 
observational volume. In the inner halo, where orbital timescales are short, (/) includes the contribution from 
multiple streams of a single progenitor. In a given local volume, we eventually will find stars belonging to 
each of these streams, while an infinitesimal volume only contains stars belonging to one individual stream. 
This is known as phase mixing. It is caused by the fact that for any individual stream, the phase-space regions 
with initially high density stretch out and get thinne r; the coarse-gra ined density decreases with time until 
a final equilibrium state is achieved (e.g. Figure 1 in lTremaine II 1999 ). Phase mixing means in practice that 
in a given volume the stellar density of tidal debris from a satellite decreases and its velocity dispersion 
increases with time, so that individual streams can no longer be separated in velocity space. The broadening 
in the velocity distribution is caused by the initial differences in the orbital phases between stars stripped at 
different epochs. In the disk for example the broadening in the velocity distribution is most prominent in the 
vertical W-component, because the vertica l frequency with respect t o the disk is shorter than the horizontal 
frequencies. For the angular y-component. lHelmi and White I (Il999h calculated that complete phase mixing 



would give rise to an observed velocity dispersion of 50-150 km s ', depending in proportion on the initial 
velocity dispersion of the satellite. 

The problem of phase mixing can be avoided, if one switches to integrals-of-motion-space (loM-space). 
An integral of motion I{r, v) is a function of the phase-space coordinates that is constant along any stellar 
trajectory in phase-space. It does not depend on time: 

I(r(t), v{t)) - const. (6) 

An integral of motion is called isolating if it isolates the points on a star's trajectory from neighboring points 
in phase space, contrary to non-isolating integr als which are infinitely many-valued (for more details see e.g. 



Ollongren Ill962t iBinnev and Tremaine II2008I) . Classical isolating integrals of motion include the energy E 
in every static potential, the three components of the total angular momentum in a spherically symmetric 
potential, or one of its components, e.g. L,, in an axisymmetric potential. According to Jeans theorem, in a 



steady state any function of n integrals of motion automatically solves the collisionless Boltzmann equation 
(|5]l, because 



><-») "-M'Zlf-"- 



dt 



(7) 



The coarse-grained phase-space density - written as a function of the integrals of motion - evolves asymp- 
totically to a state that satisfies Jean's theorem. Therefore, a subpopulation of stars initially confined to a 
small phase-space volume will remain clumped in the loM-space (Figure|2]l. 




••"■-•"■ ■.■.».;-!-3EWi- 







■1000 1000 
L, [kpc km/s] 



2000 
L [kpc km/s] 



3000 



Fig. 2 Simulated (£, |L|,L;)-distribution of stars that have been stripped from 33 different (color-coded) 
satellites and observed 12 Gyr within 2 kpc from the sun. Simulated error s as expected for the Gaia mission 
have been added to the data. Same data as in lHelmi and de Zeeuw I (120001 Fig. 4), except for the smaller size 
of the observational volume. 



Additional ins ights can be gained from studying the dispersion of stream particles in action-angle co- 
ordinates ( J, 0) jTremaine Ill999t iHelmi and White Ill999t iBinnev and Tremaine 1120081) . In Hamiltonian 
mechanics, actions and angles form a special set of canonical variables that lead to a simple form of the 
canonical equations of motion: 



d0 

dt 



dH(J) 
dJ 



= fi{J) 



dJ 

dt 



dHiJ) 
d& 



= 0, 



(8) 



where the components of fi{J) are the orbital frequencies. Neglecting the influence of the satellite's poten- 
tial on the individual stars, the actions J of a star remain constant, while the angles increase linearly with 
time. The cloud of escaped stream particles therefore spreads in only three of the six phase space dimensions 



at rates which can be shown to depend on the eigenvalues of the Hessian matrix of the Hamiltoniar|3 and the 
initial spread in the actions. For Stackel potentials that resemble disk-like flattened axisymmetric potentials, 
it can be shown analytically that the actions are functions of the three isolating integrals E, L, and 73. Again 
it follows that the stream stars remain clumped in the space of these integrals, even if the potential underwent 
slow changes in the past. The reason is that usually two or three action variables are adiabatically invariant 
(like J^ - L, in an axisymmetric potential). 

Attempts to use the action-angle variables di rectly to search for fossi l remnants of past accretion events 
have until now only been made sparsely (e.g. bv IChiba and Beers luOOOl) . The reason might be that action- 
angle variables are calculated readily only for integrable, time-independent potentials, while for real stars 
in the unknown Galactic potential they have to be c onstructed via sophisticated numerical methods (e.g. 



Binnev and Spergel ' 1984; M cGiU and Binnev 



1990l). Theoretically, the space o f orbital frequencies seems 



to be a promising one tMcMillan and Binnev luOOSt lOomez and Helmi luOlOah . In a static potential, stars 
stripped at time t — and observed at time t can be identified as distinct clumps - each representing an 
individual stream and separated by regular distances 6Qi that directly measure t through SQi — ^ (Figure|3]l. 
The reason is that any observational volume, such as the solar neighborhood, only picks out those stars that 
have moved a certain amount in an angular direction, e.g. azimuth 0, plus or minus an integer number of 
completed angular periods, like rotations about the Galactic center. The number of these integers is given by 
Qit/ljT, so the number of patches in !7-spaccl increases proportional with t^. Meanwhile, they get smaller 
proportional to r^, because the width in one dimension is given by AOj/t, where zl 0,- is the range in angle 



comp o nent ; which in turn is dete rmined by the size of the observational volume (IMcMillan and Binnev 



20081) . iGomez and Helmi I (l2010al) have shown that also for a time-evolving host potential and even a N- 



body simulation of the accretion of a satellite onto a live disk, the clumps in frequency space are conserved, 
although their regularity breaks down in the latter case. In addition, the timescale of growth of the host 
potential seems to leave its imprint as a curvature in the a lignment of the p atches - a clear advantage over 



the classical integrals E and L^, that accord ing to Penarr ubia et al. I (120061) can only constrain the present 



Galactic potential. Recentlv. iGomez and Helmi I (1201 Obi) applied the stream search in frequency space to a 
mock Gaia catalog that contained tidal debris from 42 disrupted satellites. In their simulations, these authors 
inclu ded observational er rors expected for Gaia as well as a time-varying potential of the halo, bulge and 
disk. iGomez and Hekni I found that only for small observational errors (relative parallax errors less than 
0.02), a discrimination between adjacent streams in frequency space and a clear determination of the time 
since accretion are possible. In their simulations, they were able to compute accretion times for ~ 30% of 
the identified satellites, accurate to 15-25%. The estimated times are always lower limits and the difference 
to the actual time since accretion depend on the form of the gravitational potential and its evolution. 

The actions J, besides having the advantage of being constant and adiabatically invariant once the influ- 
ence of the satellite's gravitation has ceased, also allow the identification of individual streams as different 
patches ( IMcMillan and Binnev 120081) . However, the patches are not aligned along lines of constant 7, like in 
the case of the i3,s (in a s tatic potential), because the f3,s depend on all three 7,s, so that the constraints on f2 
also affect J^. Therefore. iGomez and Helmi I (l2010ah found that individual streams are less clearly separable 
in (J,; y0)-space, comparable to the space of classical integrals (E, LS). However, the large expected number 
of streams passing the solar neighborhood given by equation (JUi could make it difficult to detect individual 
streams from a given satellite, especially when observational errors are considered. Therefore, in practice, 
it might be better to use the calssical integrals of motion (E, \L\, L.) to first detect tidal debris, because in- 



^ The elements of the Hessian matrix are given by Hij = gj.gj. - 

** Usually, the vertical frequency Q- is omitted because Q- » Qr V Q. 2> Q0, implying that the number of streams in the ^-direction 
is much larger, and so after some time, much harder to see unless one has very high numerical resolution (or the orbit is elongated in 
the z-direction). 
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Fig. 3 J7-space distribution of stream particles lost from a satellite 
without modeled self-gravity on an orbit in a spherical static poten- 
tial. Grey dots show the distribution of all satellite particles after 10 
Gyr, whereas the black dots only show the particles inside a sphere 
of 4 kpc radius. Note the alignment of the patt erns along lines of 
constan t frequencies Q, and i3^. Adopted from iGomez and Helmi I 
)2010al. Fig. 5). 



dividual streams of a satellite roughly contribute to the same clump in loM-space, raising the chances of 
identifying the individual s atellites. Afterwards, i7-space could be used to compute the time since accretion 
dGomez and Helmi 11201 Obi) . 



Indeed, the cla ssical integrals of motion have been used successfully for finding t i dal debris in the so- 
lar neighborhood dHelmi et al.lll999l: Ichiba and Beers II2OOOI: IFe Fiorentin et"ai]|20(35l: iKeplev et al 1 12007 : 



Smith et al. 112009b : however, some limitations have to be considered (pointed out bv iHelmi et al.l 120061) : 
First, the total angular momentum \L\ is only strictly conserved in a spherical potential which is not the case 
for the Milky Way. However, this seems to be a minor point, because even for stars that move in axisym- 
metric flattened potentials, \L\ is approximately conserved and simila r to the t hird integral 73, for which no 
analytic expression exists ( Binnev and Tremaine 2008, § 3.2.2). Also. lChiba and Beers (2000) showed that 
the distribution of halo stars in the space spanned by isolating integrals of motion in an aspherical Stackel- 
type potential can be closely mapped into the loM-space of a spherical potential. A second concern is that 
continuing accretion and merging events would have changed the potential steadily and altered the energy 
of an orbiting satellite. This emphazises the advantages of adia batic invariants such as the actions (among 
them J0 - L-). However, the dumpiness of the local loM-space (Morriso n et al. l2009tlKlement et al. l2009h 
suggests that (i) violent relaxation might have had a less smoothing effect on the signatures of the inner halo 
stars than expected and/or (ii) only a couple of early major mergers have contributed to the inner halo, in 
agreement with predictions from ACDM simulations and a direct effect of dynamical friction acting primar- 
ily on the most massive satellites (e.g. Helmi et al. 2003, and § l2.2l in this review). Both Dynamical friction 
and the gravitational influence of the satellite also lead to the problem that stars lost at different passages end 



up having different energy levels, in this way creating multiple signatures in the {E, L,)-space (Meza et al. 
|2005: Choi et al. 2007) . Finally, to calculate the energy, one needs the full 6D phase space information 
(r, v), which can be observed, plus the potential <P(r), which is desired, but not directly observable. 

T herefore alternative "effective" integrals of motion h ave been proposed (IHelmi et alj2006tlArifvanto and Fuchs 



2006t iDettbarn et al. 1 120071; iKlement et al. 1 120081 l2009h . The idea behind such approaches is to find pa 



rameters tha t chara cterize stellar orbits and can be easily expressed through the observables. For example, 
Helmi et al.l (l2006h proposed to look for stellar streams in a space spanned by the apocenter, pericenter 



and angular momentum L. (the APL-space). Moving groups then cluster around lines of constant eccen- 
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tricity . However, to calculate t he eccentricity, a gues s of the true potential is req uired. lArifyanto and Fuchs 
(l2006l) . lDettbarn et al. IfcoOTh . lKlement et al. I (I2OO8I) and lKTement et al. I (12009 ) followed a similai- approach 
in the sense that they tried to find stars with the same orbital eccentricity. Their strategy in findin g nearby 



stellar streams in velocity space is based on the Keplerian approximation for orbits developed by iDekker 
(119761) ; in this approximation, the potential 'P{R) is expanded with respect to ^, and the effective potential. 



0ef[(\L\,R) = <1>{R) + \L'^ l{2R^), takes on a form similar to the one in the Kepler problem. The asumption 
is that stars in the same stellar stream move on orbits tha t stay close together, which is justified by numer- 
ical simulations of satellite disruption (IHelmi et al.ll2006l) . Further assuming a spherical potential and a flat 
rotation curve, these stars should form a clump in the projection of velocity space spanned by 



V = arctan( — ), 



W 



Vaz = 



Vaz = 



V (V + Vlsr)' +"H^ 
-V(^+^lsr)2 + W^' 

Vf7^T2(VL^rTKl)2 



v< 180° 
v> 180°, 

v< 180° 
v< 180°. 



(9) 
(10) 

(11) 



V is the angle between the orbital plane (which is fixed in a spherical potential) and the direction towards the 
North Galactic Pole and ranges from 0° - 180°. Stars with v > 180° are treated as stars moving on retrograde 
orbits in a plane with inclination angle v - 180°. Vaz is related to the angular momentum \L\ and defines 
the guiding center radii of the stars, Rq - Rqjt^, in the plane labeled by v. Vae is a measure of a star's 
eccentricity e through 

e = -^ Vae (12) 

V2yLSR 

(eq. 1 1 in lKlement et al. Il2009h . Also, it can be shown that V j^ is adiabatically i nvariant, because it approx- 
imately relates to the radial action integral Jr (eq. 13 in i Klement et al. Il200 ^n Equations (l9ll- (fTTl l fully 
define an orbit in the Keplerian approximation, because its turning points are given by 



Ro 



1 +, 



and 



Rn 



Ra 



1 -, 



(13) 



Dekker I (119761) has shown that her approximation holds for eccentricities up to e = 0.5 (Vjg ~ 150 km s '), 
although in practice, the clumping in the space spanned by the quantities (l9b-(fTT]) seems to remain even for 
stronger deviations from circularity (FigureHJ. 



3.2 The use of Velocities and Positions 



As noted above, stream stars passing the solar neighborhood have very similar azimuthal frequencies or, 
equivalently, angular momenta. Therefore these stars show a narrow distribution of azimuthal velocities. In 
this respect they behave exactly like the classical moving groups in the solar neighborhood, that emerge 
either from dissolved open clusters or through dynamical resonances with the non-axisymmetric potential. 



^ In the formalism of lArifyanto and Fuchs I J2006r) and iKlement et al.1 J2008l) . slightly different notations have been used, which 
follow easily under th e ass umption of motions in the orbital plane z = or, equiyalently, v = 90° only. Note also that in the works of 



iDettbam et al. I )2007b and lKlement et al.1 12P09) only the upper definition of equation jilt has been used, which still allows identifi- 
cation of substructure, but not the calculation of eccentricities from equation {12). 
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Fig. 4 (v, Vaz, VziE)-distribution of the same stars as in Fig. [2] The same color coding has been used. The stars 
have been binned into 30° wide v-sHces, as indicated in each panel, and displayed in (Vaz, Kie). Note that 
these quantities are derived in a spherical potential, although the true potential in this case includes a disk 
(IHelmi and de Zeeuw luOO u, Model I). Also note the symmetry of multiple streams from the same satellites 
with respect to v = 90°, which follows from the symmetry in the (V, W)-distribution of phase-mixed tidal 
debris (Fig. |5]below) and equation (|9]l. 



However, moving halo groups often ha ve a distinct banana-shaped (t/, V) velocity distribution, which is 
a con sequence of their eccentric orbits JHelmi and White Ill999t iHelmi et al.ll2006t IVillalobos and Helmi 



20091) : The slightly different orbital phases of the stars in the observed volume can result in some stars 
being at their apocenters (L'^ = 0), while the others either move away from (J/ > 0) or towards (t/ < 0) it 
(Figure lU. For the same reasons, if the pericenter of some stars lies in the observational volume, a similar, 
but switched and stronger curved banana pattern results (because stars move fastest near their pericenters). 
Compared to in situ thick di sk stars, accreted stars mostly occupy the outskirts of the ({/, W)-distribution 
dVillalobos and Hekni lEooJ) . 

An even clearer separation between tidal debris and metal-poor disk stars can be obtained by observ- 
ing stars out to large cylindrical radii R. Assuming a constant rotation curve implies L, oc R for disk stars, 
whereas stream stars display a nearly constant angular momentum that depends mainly on the initial condi- 
tions of their progenitor's orbit (like the orbital inclination), but not on R. Therefore, by plotting R versus L,, 
th e disk stars discriminate them selves from the stream stars by their trend of rising L- as R increases (Fig. 6 
in lViUalobos and Helmi lE009l) . 

Finally, heliocentric line-of-sight velocities vios, which are much easier to obtain than full space motions, 
can be used to detect stellar streams against the disk star background. In this case, circular rotation of the disk 
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Stars implies a sinusoidal dependence of vios on galactic longitude ZQ while again the stream stars dominate 
the wings of the vins d istribution and should become apparent after subtraction of the mean disk rotation 
JVillalobos and Helmi ||2009.) . 



4 Moving halo groups 

4.1 Current state of research 



4.1.1 The early years 



Since lProctorl (118691) it has been known that among the nearby stars some groups exist that share the same 
streaming motion, therefore called stellar streams or moving/kinematic groups. It wa s hypothezised early that 
such moving groups result from the breakup or dispersion of open clu sters (see e. g. lEggenlll958i and refer- 
ences therein), an effect that also would apply to halo globular cl usters (IOortll965b , satelhte galaxies or o ther 
accreted building blocks of the Galaxy hke the lSearle and Zinn [ fragments. In 1959. lEggen and Sandage I dis- 
covered a group of five subdwarfs within 20 pc - among them the star RR Lyrae - with large space motions 
moving towards a common convergence point. The group was named the Groombridge 1830 group after 
the star Groombridge 1830 (HR 4550), which lags the Local Standard of Rest with V ~ -130 km s'. In 
the years that followed, Eggen extended his observations of subdwarfs in number and distance and estab- 
lished the existence of high-velocity, metal-poor moving groups (e.g. Eggen 1977, and references therein). 
Two of these, the groups named after Kapteyn's star and Arcturus, have only recently been investigated in 
great detail (see § 14. 1.3b . In the early 90s, despite some false detections of halo substructure, the evidence 
for nearby kinematic halo st reams interpreted as tidal debris increased (for a review see Majewski 19931 
Sect. 5.1).|TremaineJ (119931) computed the fractional volume containing one or more tidal streams to be 
« 1 at the Solar radius, if the halo would have been built solely through merging and accretion. 10 years 
later, Gould ( 2003) should derive from statistical arguments on a set of 4588 halo stars that there are about 
~ 400 cold streams present in the solar neighborhood, whereby no stream could contain more than 5% of the 
local stars. This number is fully consistent with the predictions of equation (|4|i and has until now not been 
proved false. More possible streams were found when Poveda et al. ( 1992) analyzed the clustering of 206 
halo stars in (E, L-, [Fe/H])-space. These authors justified their approach with the argument that "...because 
of the high incidence of chaotic orbits among halo stars. ..the usual criterium of identifying group members 
by the similarity of their V-velocities breaks down." 



4.7.2 The Helmi stream 



In 1999. lHelmi et alj (Il999i hereafter H99) published their seminal paper on the discovery of two single 
streams stemming from the same progenitor, that "...had a highly inclined orbit about the Milky Way at 
a maximum distance of ~ 16 kpc, and it probably resembled the Fornax and Sagittarius dwarf spheroidal 
galaxies." The two streams are clearly separated in the {U, y)-plot, but occupy the same region in the {L,, LjJ- 
plane(~ (1200,2200)kpckms"^)Thestreamstars areonprogradeeccentricorbits(7?peii ~ 7kpc,/?apo ~ 16 
kpc) that reach distances up to \z\m-dx ~ 13 kpc above orbelow the plane. Originallv. H99 concluded that about 
10% of the local halo sta rs would belong to this stellar stream. However, this number was questioned when 
Chiba and Beers I (120001) an alyzed 728 st ars with full space motions (accurate to ~ 20%) and [Fe/H]<-1 
which were drawn from the Beers et al. I (12000) sample of 2106 metal-poor stars without kinematic bias. 



• {U- Uq) cos b cos l+(V- Vq) cos b sin l+{W- Wq) sin h; for circular low latitude orbits this implies I'los ~ {V- Vq) cos b sin /. 



14 



Although this meant a sample size three times as large as H99' s, they found practi c ally th e same number of 
stream members. The iBeers et al. I catalog was re-examined by Re Fiorentin et alj (l2005h . but now concen- 
trating on the classical search strategy for moving groups, namely pi cking out stars with strong correlations 
in velocity space. Among the 5% of the fastest moving halo stars. iRe Fiorentin et al.l found a kinematic 
group consisting of 3 stars with similar angular momenta, but somewhat lower binding energy than the H99 
streams. They hypothezised that their moving group could consist of stars from the same progenitor of the 
H99 group, but stripped at an earlier time, before dynamical friction dragged in the satellite more closely. In 
total, based on their (L^, L i )-distributions. lRe Fiorentin et al.l c ounted 16 stars out of 410 (3.9%) as putative 



members of the H99 streams, consistent with iGould r s 120031 estimate that no single stream could contain 



more than 5% of the local halo starfl In this respect, it seems that the H99 progenitor is one of the major 
contributors to the local halo - a result that has further been supported when lKeplev et al. I (120071) analysed 



a sample of 231 nearby (d < 2. 5 kpc) metal-poor red giants, RR Lyrae and RHB stars, including again a 
subset of the iBeers et al. I (120001) catalog. These authors identified 11 or possibly 12 stars connected to the 
H99 group (~ 5%) and showed that their progenitor was likely accreted between 6 and 9 Gyr ago by compar- 
ing the observed asymmetric W-distribution of the stars (8/1 1 with negative W) to simulations. In addition, 
[a/Fe] ratios for three of the H99 stars revealed that they are very similar to other solar neighborhood halo 
stars, implying formation in a massive progenitor where enrichment occured mainly through SNe II. 

The chemical s imilarity of the s t ream s tars among each other, but also with normal inner halo stars, was 
recently shown by iRoederer et al. I (120101) through determinations of abundance s or upper Umits, respec - 
tively, for 46 elements in 12 kinematically selected member stars fr om the lists of IChiba and Beers I ( 120001) . 
Re Fiorentin et al.l ( l2005h and lKeplev et al. I (l2007h . lRoederer et al. I found that the stream stars span a broad 
range in [Fe/H] from -3.4 to -1.5, but compared to typical present-day dSph and ultra-faint dwarf galaxies, 
only small star-to-star dispersions of [X/Fe] at a given metallicity have been found. The high a abundances 
(e.g. [Ca/Fe]=H-0.4) together with the trend of the «-capture abundances (e.g. [Ba/Fe] and [Eu/Fe]) to in- 
crease with metallicity point towards a scenario in which star formation in the H99 progenitor was truncated 
before SNela or AGB stars enriched the ISM. 



Four memb ers of the H99 stream with negat ive W ve l ocities also showed up when lDettbarn et al. I ( l2007h 
re-analyzed the iBeers et al. I ([2000) data using IPekkerr s 1 19761 theory of Galactic orbits an d searching for 



stream s in the space spanned by (v, V^z, V^e) (see § 13.11 ). The same search strategy was used bv lKlement et al. 



( I2OO9I) on a sample of 22,321 nearby (d < 2 kpc) s ubdwarfs with [Fe/H 1<-0.5 from the seventh SDSS data 
release with astrophysical parameters from SEGUE ( lYannv et al. 1120091) . This sample, being fainter than the 
Beers et al. I compilation, added another 2 1 stars to the H99 stream. Referred to the total number of likely 
halo stars in the SDSS sample (~ 4400), this only accounts for 0.5% of the local halo. I should note that 
th e search strategy only picked out these 21 stars from the single stream with negative W- velocities, similar 



to 



Dettbarn et al. 's work from l2007l If one assumes the asymmetry in the W-distribution in the sample of 



Keplev et al. I (120070 as representative, another 3/11-21 ^ 6 H99 st ars with positive W v e locities could 
be exp ected, presumably too few to appear in the significance maps of iKlement et al. I (120091) . iKeplev et al. 
(120071) als o found that 67% of the ir H99 member stars lay within d = 1 kpc, a ran ge which is incompletely 



sampled in lKlement et al. I (120091 Fig. 5c). On the other hand. .Smith et al. I (120091) predict up to ~ 7% of all 



halo stars within 2.5 kpc to belong to the H99 streams. Their result is based on the fraction of 12 putative 
H99 members out of 645 SDSS subdwarfs selected from a reduced proper motion diagram, again using 
the (Lj., Lj_) criterion and correcting for any selection effects d ue to the reduced prop er r notion bias. An 
explan ation for these different results could be that the method of lDettbarn et al. I (120070 and lKlement et al. 



(120091) identifies single streams from the same progenitor in different regions of (v, Vaz, V^E)-space, requiring 



' Even more so, since it could be expected that the 16 stars from lRe Fiorentin et alj are distiibuted between three H99 streams, 
the kinematic group consisting of 3 stars plus the original two streams with positive and negative W velocities, respectively. 
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a large number of stars in any single stream for a significant detection; in contrast, these single streams all 
enhance the signal in the same region of the {L,, Lx)-plane. 



4.1.3 CO Centauri, the Kapteyn and Arcturus streams 



The iBeers et al. I ( I2OO0I) catalog, besides containing the H99 stream as the most significant feature, reveals 
further substructure when plotting a histogram of the rotational velocities of the stars like i n Figure [6l The 
peak o f stars with rotational velocities of V^oi ~ 100 km s"' was investigated in deta il by Navarro et al. 



(120041) and attributed to the Arcturus stream, first described in 1971 by Egge n (.Eggenlll971[) as a moving 
group of stars with similar kinematics as the star Arcturus. Originally, Nav arro et al. I ( 120041) suggested a 
tidal origin for the Arcturus moving group, support ed by th e well-defined sequ ence of abundanc e ratios 
of member stars taken from the then available Gratt on et al. 1 12003) catalog. Ari fyanto and Fuchs I (12006), 
who later detected the Arcturus group in a compilation of various catalogs using their (Vaz, Vae) method, 
suggested an external origin, too, based on the go odness of theoretica l 12 Gyr isochrone fits to the putative 
member stars. The same suggestion was made bv iHelmi et al.l (l2006h for a g roup of stars (their Grou p 2) 



detected in the APL-space distribution of Geneva-Copenhagen survey stars dNordstrom et al. 1 12004|) and 



possibly related to the Arcturus stream (according to its V- and [Fe/H] -distribution). 

However, recent detailed abundance studies including various a and other elements throw doubt on 



the hy pothesis that the Arcturus stream is composed of a homogeneous stellar population dWilliams et al 



20091) . Instead, the putative memberqj do not differ from the surrounding disk stars. This either indicates a 
progenitor system that had to be very massive to self-enrich to [Fe/H]=-0.6 or, more likely, a dynamical origin 
of this group. Similar to the Hercules stream, this could be either caused by the 6: 1 o uter Lindblad resonanc e 
with the bar at the solar position (B. Fuchs, unpublished) or, as suggested recently bv lMinchev et al. ( 20091). 
the rel axation of disk stars after the perturbation through a massive merger ~ 1 .9 Gyr ago. iBovev et al. 



(120091) also suggest a dynamical origin due the very small cry velocity dispersion they find for the Arcturus 
stream. 

The overdensity of stars on slight ly retrograde orb its in Figure |6](yrot ~ -50 km s ' or, equivalently, 

V 270 km s ') was first noticed by lPinescu I (120021) . most prominently when the sample was restricted to 

the metallicity range -2.0 < [Fe/H] < -1.5. As these values are consistent with those of the peculiar globular 
cluster ol> Cen , which is confined to the disk (Zmax ~ 1 kpc) and has therefore lon g been suspected to be the 
surviving nucleus of a tidally disrupted dwarf galaxy dFreeman II 1 993b . IPinescu I hypothesized tha t the peak 

at Viot 50 km s ' was caused by tid al debris from that same galaxy (see also lBrook et al. 112004 ). Further 

support for this theory came from Me za et al. I (|2005), who picked out 13 likely co Cen group mem bers by 
their L, and |W| values from the compilation of ~ 150 nearby metal-poor stars from lGratton et al. I (l2003h . 
The broad (-300 km s ' < t/ < 300 km s '), symmetric and centrally under-represented t/-distribution 
of these stars indicates that most of them currently move away and towards their apocenter, with only two 
stars having |L'^| < 50 km s"'. This agrees well with the apocenter of w Cen currently lying inside the solar 
circle at ~ 6 kpc. Apart from the two kinematic outliers, all stars additionally define a narrow track in the 
([a/Fe],[Fe/H])-plane,consistentwith the chemical composition of w Cen itself and indicative of a protracted 
phase of star formation in a self-enriching progenitor system . With ~ 7%, the oj Cen group contributes to the 
local halo comparable to the H99 stream dMeza et al. II2005I) . 



^ The I Williams et alTI 120091) selection of Arctur us group mem bers by V = -100 ± 1 km s * differs sligh tly from the definition 
of the stream in other work s: <V> = -110 km s~' teggenlll99d): (V ) = -120 km s"' JNavarro et ^112004) : (V) = -125 km s"' 
JArifvanto and Fuchs Il200q) : (V) = -105 km s"' JBovev et al. 1120091) . Note that such offsets could arise through different values for 
the LSR. 
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wCen 's kinematics also put it in the retrograde region occupied by the Kapteyn moving group. 
so by O. J. Eggen after Kaptey's star which lags the LSR with ^ 290 km s~' (Table 11 in Eggen 



named 
19961 



and references therein). Recently, 16 members of this stream listed by 'EggenJ have been observed spec- 
trosco pically to test the hypothes is that Kapteyn's group in fact originates from the same progenitor as 
wCen dWvlie-de Boer et al. 1120101) . At least 14 stars have been shown to be kinematically and chemically 
([ff/Fe] plus [Cu/Fe] for 5 stars) consistent with wCen 's stellar populations, but due to their higher ener- 
gies unlikel y to stem from the glo bular cluster itself. Instead, the possible division of the stars into three Lj. 



bands made iWvlie-de Boer et al. I hvpothesize that these stars were stripped from wCen 's parent galaxy at 
different orbital passages. The progenitor itse lf would have had a l uminosity of My ~ -11 or 2 x IQ^Lq 
and a total mass of a few IO^Mq. Interestinglv. lWvlie-de Boer et al. I estimated that such a satellite galaxy is 
expected to shed several hundred more stars into the solar neigh borhood that have yet to be detected. In this 
regard the retrograde streams found recently by various authors JDettbarn et al. luOOTt iKlement et al. Il2009l 



see also TablefTJ could add further evidence for tidal debris from wCen 's progenitor system. 



4.1.4 Further halo substructure in the solar neighborhood 



Besides the H99 and u Cen streams, which seem to dominate the fraction of nearby halo stars attributed 
to tidal debris, the recent years have added several ca ndidates for further h a lo stre ams thanks to the in- 
creasing number of large samples like SDSS/SEGUE. lArifvanto and Fuchs I (120061) noted an overdensity 
of stars with thick disk kinematics rotating approximately 20 km s"' faster than the Arcturus group, but 
with very similar metallicity and age (inferred from theoretical isochrones). They confirmed this detection 
in the G eneva-Copenhagen catalog of solar neighborhood stars. Interestingly, this catalog was also ana- 
lyzed by iHelmi et al.l (120061) . who noted an overdensity of stars in the APL space in an eccentricity range 
of 0.3 < e < 0.5, with overall kinema t ics overlappi ng with the detections of Arcturus and the new stream 
found bv lArifvanto and Fuchs I ( l2006h . iHelmi et al.l further decomposed the stars in this eccentricity range 
into three groups, based on the trend of their metallicity to decrease with increasing eccentricity in a dis- 
continous fashion. Based on the metallicity distributions, they concluded that the three groups may well be 
tidal debris from three different progenitors with initial masses of ~ 4 x 10^ M©. It should be noted that the 
{U, y)-distributions of these stars indeed show the opposite trend of what would be expected if they had a 
dynamical (resonant) origin. 

The model of iMinchev et al. I (l2009b . where the streams found bv lArifvanto and Fuchs I result from the 
relaxation of a pe rturbed disk fol lowing a major merger, is also able to expla in the location of th e moving 
group detected bv lKlement et al. I (12008.) in the first RAVE public data release dSteinmetz et al.ll200 6). which 
is at y » - 160 km s~' . The original detection of the stream in RAVE data was based on the assumption that 
the vast majority of RAVE stars are dwarfs/ subdwarfs in order to estim ate distances thr ough the photometric 
parallax method. This assumption was shown to be doubtful; in fact, ISeabroke et al. I (12008.) derived from 
reduced proper motions that ~ 44% of RAVE stars are likely K-M giants. The now available second data 
release (IZwitter et al. Il2008h . which contains spectroscopically derived \ogg estimates for ~ 21,000 stars, 
in principal allows the classification of dwarfs and (sub-)giants and a critical re-analysis of a pure dwarf 
sample (Kl ement et al. 2010, in preperation). However, evidence for the existence of the KFR08 stream was 
found by Kle ment et al. (2009) in their sample of SEGUE subdwarfs at a confidence level of ~ 99.7%. 
Furthermore. . Bob vlev et al. ■ (.2010.) identified 19 candidates of the KFR08 stream among a sample of F and 
G stars for which accurate d istance estimates (o-„/7r < 0.15) were available from the updated Hipparcos 
catalog ( van Le euwen 1 120071) : they derived a stream age of 13 Gyr from stellar isochrones and a mean 
metalliciy of ([Fe/H]) = -0.7 with standard deviation 0.3 dex from Stromgren uvbyP photometry. The 
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c hemical and chronolog ical homogeinity, as well as the banana-shaped ( U, y)-distribution of the stars (Fig. 8 
Bobvlev etlTlbOK 



in 



favor a tidal rather than a dynamic origin of the KFR08 stream. 



Ot her moving gr oup candidate s without known progenitor systems have been identified bv'D ettbarn et al. 
(l2007b and Klemen t et al. I (l2009h in the space of effective integrals of motion (v, Vaz, V ae) (§ 13. It . They are 
listed in Table 1 in th e next section. One of the features detected by iDettbarn et al. I their 'S3' overden- 



sity, was confirmed by ^ Klement et al. I and affiliated to a much larger group of stars that also show a spatial 
coherence with similarity to an overdensity detected by Ijuric et al.l ( l2008h at {R,z) ~ (9.5,0.8) kpc. The 
stream members move on slightly retrograde orbits, are currently near their apocenters and have a [Fe/H]- 
distribution characterized through (Fe/H]) ^ -1.6, cr[pe/H] = 0.4. Although this points towards a possible 
relation with w Cen , the high W velocity components and spatial concentration of the stars argue against 
this hypotheses. A detailed chemical study would be helpful in solving this issue. 



4.2 A census on the currently known halo streams 

In Table [1] I summarize all currently identified halo streams within 2 kpc from the sun, ordered by the v^ 
velocity components or, equivalently, their L. values, v^ is the rotational velocity in a cylindrical coordinate 
system and for nearby stars can be approximated by v^ x; y + Vlsr + V©. I kept the names originally given 
to the individ ual streams by the authors ci ted at the end of each row. I have omitted to list Arcturus and the 
AF06 stream 1 Arifyanto and Fuchs 2006), because their kinematics and chemistry more closely resemble 
dynamical thick disk streams (Williams et al. 2009) . 

From Table[T]it is obvious that some of the individual groups might be related in the sense that they might 
stem from the same progenitor. For example, many of the retrograde streams in the region -100km s ' < 
V0 < overlap in their rotational velocity and [Fe/H] abundances. As discussed in § 14.1.31 this region is 
thought to be populated by multiple streams fr om fa; Cen 's progenito r, and it seems very likely that many of 
these features are indeed related. In particular, Kleme nt et al. I (120091) akeady showed that the streams 'C2', 
'C3' and 'S3' are not only kinematically and chemically connected, but als o spatially coherent, similar to a 
small overdensity at (R,z) ~ (9.5,0.8) kpc described by Juric et alJ ( 120081) . Other str eams with a common 
origin might be the RHLS an d H99 grou ps (Re Fiorentin et al."2005^, R2 and KFR08 (' Klement et ^71120091) 
and the S3 streams from Dettbarn et al. 1 (2007) and Klement et al. (2009) , respectively. 

Furthermore, relations between halo streams and certain globular clusters have been suggested in several 
instances, based on their common position in the loM-space. In particular, relations have been suggested 
betwe en SKOa and the clusters NGC 5466, NGC 6934, NGC 7089 /M2 and NGC 6205/ Ml 3 (Smith et alJ 
2009r) or between the H99 stream and M15 (IRoederer et al. II2OIOI) . Also, as already discussed in § 14.1.31 
Kapteyn's stream has been linked to co Cen , and probably the globular clusters NGC 362 and NGC 6779 are 
associated, too (Dinescu 2002). 

To draw further conclusions about the origins of, and relations between, the sev eral streams one would 



need (i) more precise 3D velocity estimates (<5 km s"^ iHekni and White 1 1 1 9991) . which will ultimately 
become available for thousands of halo stars from the antici pated Gaia mission and (ii) detai led abundance 
analys es of the individual streams like the ones undertaken bv lWvhe-de Boer et al. I ( 120 lOl) and lRoederer et al. 
(l2010h . 



5 Conclusions 



In this paper, my goal was to summarize the current state of research on stellar halo streams in the solar 
neighborhood. I gave an overview on the processes that play a role in the disruption of Milky Way satellites 
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and the subsequent formation of stellar streams that eventually pass the vicinity of the sun. Disentangling 
and correlating the individual streams that contribute to the local halo thus gives valuable insights into the 
history of the Milky Way in the context of hierarchical galaxy formation. Various search strategies in com- 
plementary data sets should be used for this task in order to increase the confidence levels on the existence 



of ind i vidual streams. Further, spec tro scopic studies of the v arious streams in the style of lWilhams et al. 



( l2009b . lWvlie-de Boer et al. I ( 120 101) or lRoederer et al. I (120101) are necessary to infer clues about their origin 
(dynamical vs. accreted), properties of any progenitor system and eventual relationships between individual 
streams. To fully resolve the hundreds of expected individual halo streams, 3D velocities out to at least 2 
kpc have to be measured with accuracies better than 5 km s"'. In addition, estimating their time since ac- 
cretion requires precise measurements of parallaxes good to 1-2% dGomez and Helmi I i2010b). Ultimately, 
Gaia will deliver these informations, but - as this review hopefully demonstrated - in the meantime, many 
discoveries and conclusions are already possible with currently available data and facilities. 
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Fig. 5 Simulated velocity space distribution of satellite stars 13.5 Gyr after becoming unbound to their 
progenitor system. The figure only plots stars within a 3 kpc wide box centered on the progenitor's pericenter 
(r - 6.0 kpc, left panel), apocenter (r = 37.0 kpc, right panel) or within (middle panel). Note the banana 
shape of stars near th eir apo- and pericenters. Note that in the solar neighborhood vr ^ U, v^ ^ V and 
V, = W. Adopted from lHelmi and White I (Il999l Fig. 5). 
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Fig. 6 Histogram of rotational velocities Vrot for all stars from the 
IBeers et al. I 120001) catalog of metal-poor stars (dashed histogram). 
The distribution is decomposed into three Gaussian components repre- 
senting the thin disk, thick disk and halo as labeled. Note the excess of 
stars lagging the LSR with ~ 120 km s"' (Arcturus stream) and ~ 270 
km s"' (a)Cen stream). 



Table 1 . Current census of solar neighborhood halo streams 
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Stream 



(kms"') 



(kms"' 



([Fe/H]> o-[F,,Hi 



References 



C2 


-75 
-48" 


24 
14 


-1.6 
-1.6 


0.4 
0.4 


53 
33 


Klementetal. (2009) 


S3 


Klement et al. (2009) 


Kapteyn 


-46'' 


63 


-1.5 


0.3 


14 


Wvlie-de Boer et al. (2010) 


wCeii 


-30 








56 


Dinescu (2002) 


C3 


-24 


14 


-1.7 


0.4 


44 


Klementetal. (2009) 


CI 


-16 


9 


-1.5 


0.2 


32 


Klementetal. (2009) 


s, 


-9" 


10 


-1.9 


0.5 


10 


Dettbarn et al. (2007) 


S. 


-7 


4 


-1.9 


0.1 


4 


Dettbarn et al. (2007) 


Si 


21 


2 


-1.7 


0.5 


4 


Dettbarn et al. (200T) 


SKOa 


43 


25 


-2.0 


0.2 


6 


Smith et al. (2009^ 


R2 


59 


15 


-1.4 


0.3 


19 


Klementetal. (2009) 


KFR08 


69 


11 


-0.7 


0.3 


19 


Bobvlevetal. (2010) 


Groombridge 1830 


71' 


18 






5 


Es2en and Sandase (1^59) 


RHLS 


99 


25 


-2.0 


0.2 


3 


Re Fiorentin et al. (2005) 


H99 


140 


33 


-1.8 


0.4 


33 


Kepley et al. (2007); Klement et al. (2009) 


C4 


173 


9 


-2.3 


0.3 


20 


Klementetal. (2009) 



Note. — N is the number of total stream members given in the references, (v^) and ([Fe/H]) are the mean rotational velocity and 
metallicity of the stream, and tr,. , , o"[F;^/h] the corresponding standard deviations. 

'^From the discrepancy between the Vip velocities of the S3 stream from lDettbam et ^7ll2Q07t) and lKlement et al.1l200^ it follows 
that both authors might have described two different, but maybe related streams. 

^The V values adopted by Wvlie- de Boer et al7l l20I0) are based on modern parallaxes, proper motions and radial velocities and 
differ from the more tightly confined V values originally given in Essen 1 1996). 

lEggen and Sandagel |1959l) used Vq - 17km ; 
by later studies. 



, while I applied 5.2 km s ' . Note that this group has not yet been confirmed 
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